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We  are  interested  in  solving  the  coupled  set  of  equations 

a!’  -»Ct)y 
3^  • MtU 


For  the  problems  we  are  interested  in,  these  equations  do  not  reduce  to  any 
known  special  functions.  These  equations  arise  in  the  study  of  combat  models 
where  x and  y are  opposing  forces  and  a(t)  and  b(t)  are  the  attrition 
rate  coefficients  related  to  the  effectiveness  of  weapons.  Lanchester  [1] 
first  looked  at  these  types  of  models  for  aircraft  battles,  using  constant 
coefficients  (in  which  case  the  solution  is  in  terms  of  exponentials)  and  so 
these  types  of  equations  are  referred  to  as  Lanchester  equations.  There  are 
also  special  cases  of  the  coefficients  for  which  the  solutions  to  (1)  are 
expressible  in  terms  of  the  generalized  Airy  functions  [2]. 


We  are  interested  in  problems  where  the  coefficients  a(t)  and  b(t) 
are  positive  and  continuous.  In  that  case  we  expect  that  we  can  find  solutions 
which  are  linear  combinations  of  monotonic,  exponential  like  functions.  If 
we  have  such  a set,  then  the  analysis  of  the  forces  necessary  to  guarantee 
a win  for  one  side,  say  x , is  very  much  simplified.  So  we  are  interested 
in  determining  exponential  like  solutions  for  (1),  and  error  bounds  for  the 
solutions. 

We  convert  (1)  into  a single  second  order  differential  equations  for  x(t) 


‘ {a_1  3t}hI  " a(t)b(t)*  = 0 

x(0)  - xQ  , a'](0)  (0)  * -yQ 


(3) 


me* 


- 


B 


> 


Equation  (2)  is  still  awkward.  Let  us  define  two  useful  quantities,  the 
intensity  of  combat  I (t)  and  relative  fire  effectiveness  R(t)  , 


I(t)  = /a(tjb(t}' 
R(t)  = 


(3) 

(4) 


A change  of  independent  variables  will  simplify  (2).  Let 

t 

x = / I(s)ds  , 

then  (2)  becomes 


(5) 


d2x  D-1  dR  dx  _ n 
77  " R ’ x - 0 . 

dx 


(6) 


We  now  reduce  (6)  to  the  standard  Liouvi 11  e-Green  form  by  the  usual  change 
of  dependent  variable 


W = — . 
/R 


(7) 


So  (6)  becomes 


d2W 


l)"jw  = 0 . (8) 

Then  the  standard  estimates  for  (8)  say  that  there  are  two  solutions  to  (8) 


7? 


1 + rt-L  (1//R)| 
dx^ 


W * c-je'1  + c2eT  , 


i.e.  there  are  two  solutions  to  (2) 


x(t)  = \/ 1 {rrhe  ° 


| / I(s)ds  -/  I(s)ds  ) 


+ c2e 


1 


(9) 


What  is  the  error  in  using  (9)?  A number  of  error  estimates  are  available 
[3],  [4],  it  turns  out  that  they  are  not  quite  good  enough  for  our  purposes. 
The  estimates  depend  upon  whether  one  is  interested  in  the  dominant  solution 
or  in  the  subdominant  solution. 

(4) 


Looking  first  at  the  dominant  solution,  eT  , we  take 


w-j  = eT[l  + h i ( t)J 


Then  the  error  hj(t)  satisfies 


h^  + 2h-j  = (1  + h^ ) 


* = /R -1-  (1//R)  . (12) 

dx 

Converting  (11)  to  an  integral  equation  with  zero  initial  conditions  we  have 
hl(x)  = i/T[l  - e"2(T-s)]*(s)[l  + h^sjjds  . (13) 

Our  error  estimates  are  based  on  Will et' s generalization  of  Gronwall's  lenina. 


Lernna  1 . Let 


and  let 


h( x)  = / K(s»t)y(s)[J(s)  + h(s)]ds 
0 


I K(s, t) | s Q(s )P( x)  , 


|h(T)|  <;  P(t)  / | J(s)  |Q(s)  l’i'(s)  ( exp  j/  P(z)Q(z)  |*(z)|dz|ds  . (14) 

Proof:  See  Will et  [5]. 

An  immediate  corollary  is: 

Corollary.  Let  M(z)  = Q(z)  max  { | J(z) | ,P(z)}  then 


|h(t)|  s P(x)|expj^  M(s)|*(s}|dsj  - 1 ] . 


Proof:  Integrate  (14),  using  M(z)  . 

In  the  case  of  equation  (13) 

K(s,t)  = \ (1  - e~2(x  • s>)  s \ (1  - e*2x) 

J(s)  = 1 . 


(5) 


•/» 


% 

f 
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Thus  we  can  take  P(t)  = 1 - e and  Q(s)  = j . Then  lemma  1 gives 
Theorem  1 . The  dominant  solution  to  (8)  satisfies 


W-|  = eT[l  + h-|  (t)J 

»T  1 + (1  - e"2t)(exp^  j U(s)|dsJ  - 1) 


As  for  the  subdominant  solution,  we  take 


W2  = e“T[l  + h2(x)] 


and  then  h2  satisfies 


h,  - 2h,  = +(1  + hj 


which  converts  to 


h2(x)  = f I[e2(T"s)-  l]  4s)[l  + h2(s)]ds 


Now  we  have 


K(s,t)  = i e2(x‘s)-  1 s i (e2x  - 1 )e"2s 


Then  we  have 


Theorem  2.  The  subdominant  solution  to  (8)  satisfies 

W2  = e"Tt1  + h2^T^ 

s e‘T  1 + (e2t  - 1)  / e'2s  exp  / (1  - e_2z)  j^i^-ldz  ds  . (19) 

0 s -I  QED 

Observe  that  we  did  not  use  the  corollary  on  (19). 

We  can  see  that  our  estimates  depend  very  strongly  on  ip  , given  by  (12). 
For  the  Lanchester  equations  (1),  the  coefficients  a(t)  and  b(t)  , for 
many  applications,  can  be  expressed  by 

a(t)  = k (t  + c)u 

a (20) 
b(t)  = kfa(t  + c + a)u 


(6) 


where  c is  the  "starting"  parameter,  and  A is  an  "off-set"  parameter 
(so  y can  have  a different  firing  range  than  x ).  These  are  referred  to 
as  the  "power  rate"  coefficients.  For  future  use,  let 


B * /k  k, 
a b 


r _ U + U + 2 

o - * 


We  first  consider  the  case  of  no  offset,  i.e.  A = 0 . Then  (see  (5)) 


= - 


u - u (u  + 3u  + 4] 
16  (fir  + C6)2 


We  see  that  for  -1  < u < u that  tp  is  negative,  and 


j^exp  / | ds j - 1 

■ -1  * ■ 


r -=  W1  + ^ + 4> 

and  D<j  * exp  (y/A  ) . 


Then  (16)  becomes 


Wj  s eT(  D-j  exp  (-y/(A6  + St) ) + e - D-j e 3 1 exp  (-y/(Aa  + 5t)) 


as  an  error  term  for  the  dominant  solution. 

For  the  subdominant  solution  the  corollary  gives  too  crude  an  estimate. 
Even  with  using  the  lemma,  the  error  estimate  is  not  good.  Integrating  by 
parts  we  get 


(7) 


*2  = e"  T [1  + h2(  x)J 

s e*  1 [1  + (e2T-  I)(l  - exp[7  0 - e-2z)|f|dz] 
+ ^|f|exp[7'(l  - e-2z)|ti|l|<JZ]dsy 


Working  this  out  we  get 


w2  < e 


.^+0(exp(_x_))] 


This  is  not  a very  exciting  bound,  but  we  are  unable  to  do  any  better,  using 
the  Will et  result. 

We  can,  however,  alter  the  Wi 11 et  result.  The  subdominant  solution  error 
also  satisfies  the  equation 

h = I /“[l  - e"2(s  ' t)>(s)[1  + h2(s)Jds  . (26) 


We  now  prove 

Lemma  2.  If  h(  t) 

satisfies 

h(-u)  s / K(s,  t)ip(s ) [J(s ) + h(s)]ds 

T 

(27) 

and 

L 

| K(s,  t)  | £ P(x)Q(s) 

(28) 

where  K,  i|i,  and 

J are  all  s 0 , then 

h(-r)  s; 

P(t)  / ♦(s)Q(s)J(s)exp|7  P(a)Q(a)ip(a)dcjds 

(29) 

Proof.  Eq.  (27)  can  be  written,  using  (28) 

<;  / Q(s)*(s)J(s)ds  + / Q(s)<|>(s)h(s)ds  . 

r'T'  T I 


Di  fferentiating 


I 

(~)  <;  -Q(t)^/(t)J(t)  - Q(T)*(T)h(i) 


Therefore 


(£)  + Q*P(£)  * -Q*J 


Integrating  we  get 


s / Q(sH(s)J(s)exp^/  P(<j)Q(a)i|>(a)dajds 
+ K exp^/  P(a)Q(o)>Ka)daj  . 


But  h(«)  = 0 , so  K = 0 . 

For  equation  (26)  we  have 


K(s.x)  = J-[l  - e-2(s  ' t)]  < 1 • [ 

h2  * / exp  j^/  | (a)dajds 

= exp  I / da  | - 1 . 


This  is  the  same  bound  found  by  Olver  [3], 

For  the  case  of  the  power  attrition  coefficients  with  no  offset  we  have 


w2  s e 


*exp(-T^_\ 
Vr  + 5t/ 


which  is  a fairly  reasonable  bound. 

For  the  case  where  y = u = 1 , the  linear  case,  we  can  also  get  some 
explicit  results  with  offset.  That  is 

a(t)  = k (t  + c) 

d 

b ( t ) = kb(t  + c + A)  . 

Then 

t,3A£pT7-tn(n*^)j  I 


(9) 


where 


Then 


n - i + 2(t  + cj 
n - i + a • 


■k  } ^(c)da 


T0 


1 U(t  + c)3  - 6A(t  + c)2  - 1 2A2 ( t + c)  - 7 A 


1 6BA 


[(t  + c)(t  + c + A)J 

A2c 

37? 


3/2 


4c3  - 6 Ac2  - 12A2c  - 7 A3 


£C(C  + A)]' 


while 


j / v(a)da 


1 


16BA 


_ 4(t  + c)3  - 6A(t  + c)2  - 12A2(t  + c)  - 7A3 
[(t  + c)  ( t + c + A)]3/2 


Both  integrals  are  monotone  increasing. 
For  the  general  case 


j ! <^(s)ds 


-! 


n ( (u  + |v2)(n3  + An  + A2/4)  - (u  + u2/4)(n^  - An 


2,/ka'<b  n0 


(n  - A/2)2  + u/2  (n  + A/2)2  + u/2 


(n  - A/2) 


1 

j uu 

1 + p/2 


(n  + A/2) 


r '+  u72 


dn 


(33) 


(34) 


+ A2/ 4) 


(35) 


where  n = t + c + A/2  . 
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